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^ ■ Abstract 

O ! It is known that the radial equation of the massless fields with spin around Kerr black 



holes cannot be solved by special functions. Recently, the analytic solution was obtained 
by use of the expansion in terms of the special functions and various astrophysical ap- 



^ . plication have been discussed. It was pointed out that the coefficients of the expansion 

I by the confiuent hypergeometric functions are identical to those of the expansion by the 

hypergeometric functions. We explain the reason of this fact by using the integral equa- 
tions of the radial equation. It is shown that the kernel of the equation can be written by 
the product of confiuent hypergeometric functions. The integral equaton transforms the 
expansion in terms of the confiuent hypergeometric functions to that of the hypergeo- 
metric functions and vice versa, which explains the reason why the expansion coefficients 
are universal. 
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1 Introduction 



One of the remarkable feature of the Kerr black hole, which is known to be unique 
black hole solution with no electro-magnetic charge|jl|, |[, is that the gauge invariant per- 
turbation of the spin 00, the 1/20 fields, the electromagnetic fields^], the gravitinoP| 
and the gravity^ can be treated by separable equations called Teukolsky equations]^]. 
It should be noted that Whiting||^ provided an analytic proof of the stability of the 
Kerr black holes by using these equations. Unfortunately, Teukolsky equation cannot 
be solved by any popular special functions since the equation has two regular singular 
points at horizons and one irregular singular point at infinity. In order to determine the 
physical quantities such as the absorption probabilities of black holes, we usually had 
to use numerical integration of the equation. However, for the purpose of investigat- 
ing the property of the scattering data, it seems useful to provide analytic solutions as 
has been achieved for Coulomb scattering. These may be used as starting points of the 
perturbation for the quantum corrections. One of physical application is the absorption 
coefficients of the black holes[0, H. Recently, the analytic methods have been getting a 
powerful tool even in the field of the gravitational wave astrophysics. |10|, [11], 



The technique for analyzing this type of equation is quite old because a equation 
similar to Teukolsky equation appeared when we consider the wave function in terms 



of spheroidal coordinates, which is known to be spheroidal wave equation[|T^, 14]. The 
solution of the equation can be expanded either by the hypergeometric functions or by 
the confiuent hypergeometric functions. Therefore it is natural that similar expansions 
has been considered even for Teukolsky equation]^, 0, ^. 

Regarding the structure of the expansion coefficients of Teukolsky equation, an inter- 
esting observation was given in Refs. 0, |^ who showed that the expansion coefficients in 
terms of the hypergeometric functions are identical to those by the confluent hypergeo- 
metric functions up to re-definition of the coefficients. This universality of the expansion 
coefficients appeared even in the case of the spheroidal wave functions ||13|, The wave 
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function of three dimensional flat space is separable in the spheroidal coordinates. The 
radial equation and the angular equation can be identified by some change of variables, 
which can be identified as a single spheroidal wave equation. The angular equation is usu- 
ally expanded by Legendre functions whereas the radial equation is expanded by Bessel 
functions. It turns out that the expansion coefficients are identical up to some normal- 
ization factor. The reason was explained by use of the fact that the angle equation and 
the radial equation of spheroidal wave equation are connected by an integral equation, 
because of which the coefficients are identical up to normalization factor [^. This fact 
indicates that there might be some transformation which connects these expansions even 
for Teukolsky equation. However we cannot obtain the integral equation by a simple 
generalization of the spheroidal functions because the construction of the kernel in the 
case of spheroidal functions crusially depends on the fact that the original space is just a 
three dimensional fiat space. In the case of Teukolsky equation, we cannot construct any 
simple system even if we associate any angle variables. Therefore, the above construction 
cannot be applied. Instead of considering an analog of such spheroidal wave functions, 
we treat the equation as an analog of Heun's equation |jl6|. Heun's equation[0 has four 
defenite singular points in the equation and Teukolsky equation can be considered as a 
confluent limit of Heun's equationjl^. We will use a principle used for the construction 
of the integral equation of the Heun's equation|T^. By constructing the integral equation 
for the Teukolsky equation, which can be regarded as a confluent analog of the kernel of 
the Heun's' equation, we will show that the expansion coefficients are universal. 

In the next section, we review the analytic expansion of the Teukolsky equation. 
There are basically three type of expansions which cannot be connected by the analytic 
continuation. 

In section 3, we will construct a integral equation of the Teukolsky equation. It 
will be shown that the integral kernel can be written by a product of the confluent 
hypergeometric functions. 

In section 4, it will be shown that the these expansions can be connected by the 
integral transformation, which implies that the expansion coefficients are universal. 
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Section 5 is devoted to some discussions. 



2 Analytic expansions of Teukolsky equation 

In the Boyer-Lindquist coordinates and in units such that c = G = 1, the Kerr metric 
is written as 

2Mr\ f AMarsin^ e\ , ,, /S 



ds^ = f 1 ^ j dt^ + I J dtd(l) - I ^ 1 



^,.2 2^/2 2 2Ma'^rsm^e\ , 2 . ^ 
-Ede^ - sin^ e Ir^ + + \d^^, (2.1) 



where M is the mass of the black hole, aM its angular momentum, E = r + a cos 9, 
and A = + — 2Mr. Klein-Gordon equation for massless fields ip in a Kerr black hole 
background can be separated by setting ip = e~*'^*e*™'^S'["(6')i?(r), and radial equation is 
given by 

, df^^^,dR\ fK^-2is(r-M)K \ , , 

A-'— A^+^— + )■ + 4isur -\]R = 0, 2.2 

dr \ dr J \ A J 

where s is a parameter called spin weight of the field, and X = E — 2amuj + o^uo"^ — s(s + 
1), and E is eigenvalue of spheroidal harmonics S'^{9). This equation has two regular 
singularities at r = r± = M ± VM2 - a2 = M±p and an irregular singularity at r = cxd. 
Setting the new variable z = (r+ — r)/2p, radial equation becomes^ 

z'iz - If^ + {s + l)z{z - l){2z - 1)^ 

"7^2 isK(2z — 1) 

+ — r + ^ - z(z - l)(8isujpz - Aisujr+ + X) i? = 0,(2.3) 

4^"^ 2p J 

where s takes integer or half integer value and K = {Ap'^z'^ — Apr^z + 2Mr+)a; — am. 
The analytic solution of this equation is obtained not by any special functions but by 
expansions in terms of special functions such as the hypergeometric functions and the 
confluent hypergeometric functions. 

If we use the expansion in terms of the hypergeometric functions, we have two in- 
dependent solutions of equation ( [^.31 ). Introducing the slightly modified hypergeometric 



function P(a, b; c; z) 



P(a,5;c;^) = lML^F(a,6;c;z), (2.4) 
Tic) 



where F{a, b; c; z) is the hypergeometric function, we can express two solutions around 
z = in the forni]^ 

oo 

R'^ = e''^\-zyil-zY J2 h':P{a,M;c;z), 

k=—oo 

oo 

R- = e^'^i-zY+'-'^il - zY htP{ak - c + 1, 6^ - c + 1; 2 - c; z), (2.5) 

fc=— oo 

where 

p = —s — iujM — IT, 5 = iujM — IT, a = 2puJ, (2-6) 
and r = {2ujM'^ — am)/2p, and 

aj^=p + 5 + s-k-v, bk = p + 5 + s + k + v + l, c = 2p + s + 1. (2.7) 

By inserting the expression (2.5) into (2.3), we find that the expansion coefficients h]ysjid 
hl'^ should satisfy the same recursion relations; 

akhl+, + PkK + 7kK_, = 0, (2.8) 

ia{s - 2iujM + k + u + l)(-2ir + k + u + l){s + 2iujM + k + v + 1) 



{2k + 2v + ?,){k + V + I] 



9 9 9 9 X., X 2aT{2ujM + is){2uoM -is) , , 

i5k= 8u^M^ - aW -E + {k + iy){k + iy+l) + \^ , . ^ , . ^, 2.9 

[k + u){k + z/ + 1) 

_ ia{s - 2iuM -k- z/)(-2ir -k- v){s + 2iujM - k - v) 

{2k + 2v -l){k + v) ■ 

The parameter v can be obtained by the convergence of the recursion relation 0, ^ or by 

using the post-newtonian principle 

Note that the solutions around 2 = 1, which represents the inner horizon, can be 

obtained by the analytic continuation of the hypergeometric functions. Two independent 

solutions around z = 1 are 

oo iZu 

Rl = e"'\-zr{l-zY FT h 7-,P{akM;ak + bk-c+l-l-z), 

fc_^r(c-afc)r(c-6fc) 



Rl = e''''{-zy{l~zY (2.10) 
X V r=- —P{c-ak,c-bk;c-ak-bk + l;l-z). 

Since these expansions are expressed as linear combinations of Ri and R2 due to the 
analytic continuation, the recursion relations for h^^ and hf" are identical to those of h^", 
and hf'. 

Other type of solutions are given by the expansion in terms of the confluent hyperge- 
ometric functions. Introducing the new variable p = —2ujpz, one type of the expansions 
is given by 0, § 

00 

R'^ = p\p + 2copr'~^ Y: gl''Gk+Av,p), 

fc=— 00 

00 

R'^ = p\p + 2u;pr'~^ Y: ^7rG_fc_,_i(r^,p), (2.11) 

fc=— 00 

where rj = —{2ujM + is) and 5 = iujM — it, and Gi{r],p) satisfies the Coulomb type 
equation 

Pg^Giiv, p) + 2^GKr], p)+\p-2r^- j G,{r^, p) = 0, (2.12) 

which is expressed in terms of Kummer's confluent hypergeometric function $(a; b; p) as 

Gi{ri, p) = ^!i^^(0^f''^^'P)'^^^ + 1 - zry; 2/ + 2; 2^p). (2.13) 

By inserting (|2.11| ) into the original equation, we find that the expansion coefficients g\'^ 
and g'^'^ should satisfy the following recursion relations; 

+ « + 7^^?Li = 0, (2.14) 



ik= 7fc, 



(2.15) 



which is identical to (|2lS 
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Moreover, if we set x = 2ujp{l — z), which is the expansion in terms of the inner 



horizon, another expansion can be obtained in the form ^5 

oo Sf 

R^, = X'ix-2u;p)-^-^ Y: ^, j-^G,^,{r^,x), 

fc_^r(c-afc)r(c-6fe) 

oo Au 

K = x'(x - 2.p)-< r(c-a!)r(c-M ^---''"- '^)- P'^^' 

where e = —s — iuoM — it. It turns out that the recursion relation for and are again 
identical to those of g\'^ and gl'^ and /i^. In this case, we cannot explain the reason of 
having identical recursion relations simply from the analytic continuation on the contrary 
to the case of hypergeometric functions. 

In the above construction, the coefficients g'^ and /i^ are identical up to normalization 
factor as 

9l - K. (2.17) 

Since /i^, g^ are expansion coefficients in terms of different kinds of special functions, 
such as the hypergeometric function and the confluent hypergeometric function, which 
have different analytic properties and different asymptotic behaviors, it is not trivial why 
the relation (|2.17| ) holds. Moreover, the identification of the coefficients g]^ and g^'' is 
also the subject which cannot be explained by any analytic continuation of the special 
functions. In the next section, we are going to explain the reason by use of the integral 
equations. 



3 Integral equation 

We are going to construct the integral equation in terms of R{z) which satisfies the 
equation 



M,R{z) (3.1) 
d'^R{z) fl 1 \dR{z)\ 



z z — Ij dz 
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+ 



isK{2z - I] 



Sisupz + Aisur^ — X 



R{z) = 0. (3.2) 



Ap'^z{z - 1) 2pz{z - 1) 
The integral transformation which maps one solution to other solution is given by 



R'{x) = J^y'{y-irK{x,y)R{y)dy, 



(3.3) 

where the function R'{x) is also a solution of equation ( p.3|) if the kernel satisfy the 
condition 

iM,-My)K{x,y)=0, (3.4) 
and the surface term of the integral 



vanishes at the end of C 17 . To find the kernel, we set new variables as 



(3.5) 



^ = xy, C = {x-l){y-l] 



(3.6) 



Then the equation (|3.4| ) becomes 



^^ + (^ + 1^ 



dK 



+ <! Ap'^uj'^^ - SpMuj'^ - Aisujp + 



{Mu + t){Muj + t -is) 



K (3.7) 



C,—--{s^\)—~Up'^X^ 



2 2, , {Mu-T){Mu-r-%s) ' 



K = 0. 



Therefore we can separate the variable, and we put K{x,y) = P{C,)Q{() so that P{i) 
and Q{() satisfy the equations 



e^ + (. + i)"^(^) 



(3.8) 



f , (Muj-T)(Muj-T-is) 

+ UpWc + ^ '-^ 



Q(C) = Ag(C), (3.9) 
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where A is a separation constant. Each equation has two independent solutions and these 
are expressed by using Kummer's confuent hypergeometric functions as 

P2{0 = e-^r''-^^ [-^-^^f^ - 1 - 2p - -2^a^^ , (3.10) 

Qi(C) = e-^C^$ [^^^{^ - 2!^' 25 + +. + 1; -2zac) , 

Q2(C) = e^'-^C'-'^ [- ^^^2~^ ~ 2^' ^ " " "2'^^) ' ^^-^^^ 
where p = —s — iMuj — it, 5 = iMuj — ir, a = 2puj and E = SMuj'^p + Aisup. We can 
thus solve the kernel equation in a general form. Namely, the kernel K{x, y) is expressed 
as the simple product of these solutions, or as the product of some linear combinations of 
these solutions such as Whittaker's function. Since A is a arbitrary constant, we choose 
A in such a way that the kernel becomes a simple form in order to evaluate the integral 
transformation easily. The integral equation maps a solution of Teukolsky equation to 
another solution. Therefore, we can transform a analytic expansion to other expansions 
by various choice of the kernel. In the next section, we are going to perform the integral 
transformation of such various expansions of equation (|2.3| ). 



4 Integral transformations 



First of all, we consider the transformation from the hypergeometric expansions (|2T 
to Coulomb expansions ( p^.llD around 2; = 0; we take Ri{y) in ( [^.5D as R{y) 



R\ = e'-\-yY{l-y f P{ak,h; c;y). (4.1) 

fc=— 00 

A suitable choice of the kernel for this transformation is as follows. We choose P{^) in 
the kernel by using Whittaker's function Mi^^^{z) as 

P{^)=r^^^^y^M,,,{-2zaO, (4.2) 



where k = {E + X)/2ia and fi = p + s/2. In order to eliminate the factor (1 — y) in the 
integral equation, we take (52(C) ( P-llj ) for (5(C)- Fo^' the convergence of the integral, 
we take A = ia{—26 — s + 1) so that Q{C) becomes 

QiO = e'-^C'"- (4.3) 
Then the integral equation is given by 

R[{x) = {const.)x-^'+^^^^{x - ly^-^e-'"^^-^^ 
» 00 
X y e^-^i/'(^+i)/2+''+^M«,^(-2mxy) hl^P{a,M:c;y)dy. (4.4) 

By taking the region of the integration as the interval from to +00 and using the 
formulae in |T^, we can evaluate the integral as 

00 

E'M = {const.)p\p + 2up)-'-' Y: hi'' 

fe=— 00 

X {Cir(2A; + 2z/ + l)T{-2k - 2iy)G-k-u-i{v, p) 

+C2T{-2k -2v- \)Y{2k + 2z/ + 2)G'fc+,(r/, p)}, (4.5) 

where Ci are constants which are independent of k. Note that if we start with i?2; 
we obtain the same result as the equation ( [4 .51) after the integral transformation by using 
an appropriate kernel. 

Note that the solution R'i{p) consists of two independent solutions which are expanded 
in terms of Coulomb type functions and they are identical to the solution ( |2.11|) . Thus 
two kinds of expansions are connected by the integral transformation. Recognizing the 
coefficients of G^^kijli P)i G-k-v-iijIi P) Qk^i Qk" respectively, the relation between h]y 
and gl"", gl" is 

gl'' ~ T{2k + 2v + l)V{-2k - 2u)hl^ ~ (4.6) 
^2. ^ r(_2fc _2u- l)T{2k + 2v + 2)hl!' ~ /^^^ (4.7) 

Therefore the fact that coefficients h^^^, g^ of different kinds of expansions satisfy the same 
recursion relations can be understood quite naturally. 
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Let us consider the inverse transformation. In this case, we start with a shghtly 
modified form of -Ri(p) in (|2.11| ) as R{y) 



k=—oo 

fc=-oo e 2 r(2A; + 2z/ + 2)p 

Our choice of the kernel is as follows. In order to eliminate the factor {y — 1), we take 
Qi{C) with A = -ia{26 + s + 1) so that 

QiO = e-^'^^C'. (4.9) 

For the convergence of the integral, we combine Pi(^),P2(0 i^^o ^^e form 

p(^)=^-(.+i)/2^_^^^(2^aO, (4.10) 

where W^f^^^^z) is Whittaker's function. Then the integral equation becomes 

R[{x) = {const.)x-^'+^^/^ix - l)^e*-(--^) 

By taking the region of the integration as the interval to infinity and using the formulae 
in ||20[, we can evaluate the integral. The solution is 

00 

^^{x) = {const.)xf{l-xYe''''= ^ qI"" {CiP{ak,hk]C]x) 

k=—oo 

+ C2 x^-'Piak -c+l,bk-c+l;2-c;x)} . (4.12) 

Thus we could perform the inverse transformaion, which completes the relation between 
two expansion around z = 0. The relation between expansion coefficients are 

hi-" - hf ~ gl^ ~ gl\ (4.13) 
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We next consider the relation between expansions around z = 1. As before, we use 
R'^ix) in (|2.16 ) as R{y) in the modified form 



oo 3i/ 

Rsix) =X'{X- 2ujp)-^-^ Yl w TfT. T-rGk+Av, X), 

,^^^T{c-ak)T{c-bk) 

= x'-\x-2iupr-' (4.14) 

r(c - ak)T{c - bk)e 2 T{2k + 2z/ + 2) 

where x = 2u;p(l — y) and e = — s — zcijM — zr. In this case, the role of y and 1 — y are 
interchanged in the integral equation. In order to eliminate the factor y, we take P{^) as 

P(0 = Ce"-^, (4.15) 

where we set A = ia{2e + s + 1) — E. For the convergence of the integral, we take Q{C) as 

Q{C) = C'^^+'y'W.,,,,{-2zaO, (4.16) 

where n' = X/2ia, fi' = 6 + s/2. Then the integral transformation becomes 

R'^' = {const.)x'{l - x)-(^+^)/2e'"" 

X /(l-|/)^-(^+i'/2+^-ie'"^(^-^)iy«,„,(-2i(T(l-x)(l-i/)) (4.17) 
Jc 

f^^^T[c- ak)T[c-hk)T{2k + 2v + 2) 

By considering the integral region as the interval to infinity, and by using the formulae 
pO| , we obtain the solution of the equation as 



m 



00 

i?r = (con.t.)x^(l-x)V'^^ E 9f 

k=—oo 



r(c - afe)r(c - bk) 



P{ak,bk;ak + bk - c+l;l - x) (4.18) 



(72(1 - x)'^"''''"^'= 1 

+ —, r=r, rT^(c - ak,c-bk;c- ak-bk + l;l- x)} . 

T{c- ak)T{c-bk) J 

Note that the same result holds if we start with R^x). 

R'^ consists of two independent solution which are expanded in terms of hypergeo- 
metric functions around z = 1. Thus two kinds of expansion around z = 1 are connected 
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by the integral transformation. The relation between expansion coefficients are 

gf - gt ~ hf ~ ht. (4.19) 



We thus find that all the coefficients are connected by the analytic continuation of the 
wave function and the integral transformation which maps the solution of Teukolsky 
equation to other solutions. 



5 Conclusion 



We have constructed the integral equation in terms of equation (|2.3|) , and solve the 
kernel in the general form as the product of special functions. By various choice of the 
kernel, we have performed the integral transformaions which connect various expansions 
of equation (|2.3|) . In all regions of z, we can relate these expansions one another by the 
integral transformation, and by the analytic continuation of hypergeometric functions. As 
a consequence, the coefficients g^^ of Coulomb type expansions satisfy the same recusion 
relations of the coefficients h^^ of the hypergeometric expansions, and they are identical 
up to normalization factors. 

In any case of integral transformations, the kernel is expressed as the product of the 
confiuent hypergeometric functions. It is this property that makes us possible to perform 
the integral transformation. 

Let us consider the case of the spheroidal equation. The equation ( p.4|) is just the 
wave equation in the spheroidal coordinates where the variable x is the radial coordinate 



and y is the coordinate representing the angle pTSf . In other words, x and y can be treated 
as dual coordinates. It was quite easy to obtain the integral equation of spheroidal wave 
functions because we know other type of separable coordinates in fiat space. On the 
other hand, in the case of Teukolsky equation, we have dealt with dual coordinates of 
the radial function which are not the spheroidal coordinates in Kerr geometry for the 
construction of integral equation. The existence of other separable variable ( p.6|) seems 
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to show some kind of symmetry of the space-time in a more wider geometry. Further 
study on this point of view seems interesting. 
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